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ABSTRACT 

The  transport  porperties  of  a  dilute  gas  are  determined  by  binary 
collision  integrals.  In  order  to  predict  the  first  density  corrections  to 
the  transport  properties,  it  is  necessary  to  consider  collision  integrals 
that  account  for  the  effects  of  collisions  among  three  molecules.  In 
this  technical  report  v/e  derive  and  formulate  such  three -particle  col¬ 
lision  integrals  for  the  coefficients  of  thermal  conductivity,  viscosity 
and  self-diffusion  of  a  gas  of  hard  spherical  molecules.  An  evaluation 
of  these  three-particle  collision  integrals  will  be  presented  in  a  subse¬ 
quent  technical  report. 
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I .  Introduction 

The  theory  for  the  transport  properties  of  dilute  gases 
consisting  of  spherically  symmetric  molecules  is  well  estab¬ 
lished.  This  theory  is  based  on  the  Boltzmann  equation, 
which  can  be  solved  using  a  procedure  introduced  by  Chapman 
and  Enskog.  The  transport  properties  are  thereby  expressed 
in  terms  of  binary  collision  integrals.  Calculation  of 
these  binary  collision  integrals  has  become  a  routine  pro¬ 
cedure,  and  tabulated  values  are  available  for  many  forms 
of  the  intermolecular  potential  [1,2]. 

The  validity  of  the  Boltzmann  equation  is  restricted 
to  the  low  density  limit.  In  order  to  evaluate  the  trans¬ 
port  properties  of  a  moderately  dense  gas  we  consider  an 
expansion  in  terms  of  the  density.  Originally,  it  was 
envisaged  that  the  transport  properties  could  be  represented 
by  a  power  series  in  the  density,  in  obvious  analogy  to  the 
virial  expansion  for  the  compressibility  factor.  Subsequent 
developments,  to  some  extent  discussed  in  AEDC-TR-69-68  [3], 
have  revealed  that  the  density  dependence  of  the  transport 

properties  is  more  complicated,  and  that  terms  logarithmic 
■  + 

in  the  density  appear  also.  That  is,  the  density  expan¬ 

sions  for  the  thermal  conductivity  A,  the  shear  viscosity  n 
and  self-diffusion  coefficient  D  are  of  the  form 

f 

For  a  bibliography  on  the  subject  see  Ref.  [4] . 
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2 

X  =  Xe+X^n+X2n  log  n+*** 

2 

n  =  Ti0+n^n+ri2n  log  n+«*«  (1.1) 

2 

nD  =  D0+D1n+D2n  log  n+*** 

where  n  is  the  number  density.  The  lowest  order  terms  X6, 
ri0  and  D0  are  given  by  the  Chapman -Enskog  theory-  In  this 
report  we  focus  our  attention  on  the  coefficients  X^, 
and  of  the  contributions  linear  in  the  density.  Empiri¬ 
cal  values  for  the  coefficients  X.^  and  of  a  number  of  gases 
were  reported  in  AEDC-TR-69-68  and  AEDC-TR-71-190  [3,5,6]. 

In  order  to  calculate  these  coefficients  from  the  theory, 
it  is  necessary  to  evaluate  collision  integrals  that  are 
related  to  sequences  of  collisions  among  three  molecules  [3] . 

Prior  to  the  systematic  developments ,  Enskog  made  an 
intuitive  attempt  to  generalize  the  Boltzmann  equation  to 
a  dense  gas  of  hard  spheres  [7] .  The  Boltzmann  equation 

is  akin  to  the  perfect  gas  law  in  that  it  neglects 
any  correlations  in  the  positions  of  the  molecules; 
furthermore,  the  velocities  of  two  molecules  which  are  about 
to  collide  are  also  assumed  to  be  uncorrelated  (assumption  of 
molecular  chaos) .  In  the  theory  of  Enskog  for  dense  gases 
an  estimate  was  made  to  account  for  the  effect  of  correlations 
in  configuration  space;  they  were  assumed  to  be  the  same  as 
those  for  a  gas  in  equilibrium.  The  assumption  of  molecular 
chaos  was  retained  for  the  probability  distribution  in 
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velocity  space. 

In  a  previous  report  estimates  for  the  complete  triple 
collision  integrals  were  presented  for  a  gas  of  hard  spheri¬ 
cal  molecules  [3,8].  It  turned  out  that  their  contribution 
differed  by  less  than  ten  percent  from  the  value  predicted 
by  the  theory  of  Enskog.  This  conclusion  was  recently 
confirmed  by  Henline  and  Condiff  [9] . 

In  this  report  we  shall  show  that  for  a  gas  of  hard 
spheres  the  theory  of  Enskog  is  recovered  as  the  first 
approximation,  when  the  collision  integrals  are  expanded 
in  terms  of  the  number  of  successive  binary  collisions 
between  the  molecules.  For  this  purpose  it. is  important  to 

V 

make  a  distinction  between  the  effects  of  statistical  versus 
dynamical  correlations.  In  the  equilibrium  case  the  first 
density  correction  to  the  radial  distribution  function  of 
particles  1  and  2  is  determined  by  those  configurations 
for  which  a  third  particle  is  overlapping  with  both  particles 
1  and  2.  That  is,  the  configurations  which  determine  the 
excluded  volume  are  those  for  which  the  center  of  the  third 
particle  would  lie  inside  the  action  spheres  of  both  particles 
1  and  2.  If  in  the  non- equilibrium  case  the  integration  in 
the  triple  collision  integral  is  restricted  to  the  same 
configurations,  the  triple  collision  integral  yields  for 
the  first  density  correction  the  equilibrium  radial  distri¬ 
bution  function  multiplied  with  the  Boltzmann  collision 


3 


AEDC-TR-72-142 


integral  [3] .  This  is  precisely  the  contribution  given 
by  the  Enskog  theory  and  it  contains,  effectively,  the 
dynamics  of  only  one  binary  collision  among  the  three  parti¬ 
cles.  Deviations  from  molecular  chaos  can  then  be  incorpor¬ 
ated  in  successive  higher  approximations  by  considering 
correlations  brought  about  by  two,  three  and  four  successive 
collisions  among  the  particles.  The  details  of  this  expan¬ 
sion  will  become  more  clear  in  Sections  IV  and  V. 

The  purposes  of  this  report  are: 

a)  to  present  a  simple,  self-contained  derivation  of  the 
three-particle  collision  integrals  that  bypasses  the 
somewhat  intricate  surface  integration  method  proposed 
by  Green  [10]  and  further  developed  in  previous 
technical  reports  [3,11]. 

b)  to  elucidate  how  an  analysis  of  the  statistical  and 
dynamical  correlations,  caused  by  the  interactions  of 
three  molecules,  leads  to  a  decomposition  of  the  collision 
integrals  associated  with  the  number  of  successive 
collisions  among  three  particles. 

c)  to  express  the  three-particle  collision  integrals  in  a 
form  suitable  for  a  quantitative  evaluation  of  these 
integrals. 

The  treatment  is  restricted  to  a  gas  of  hard  spherical 
molecules.  In  this  case  the  molecular  dynamics  reduces  to 
sequences  of  successive  binary  collisions  among  the  molecules. 
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The  analysis  is  simplified  by  the  fact  that  we  do  not  need 
to  consider  more  than  four  collisions  as  demonstrated  by 
Hoegy  and  Sengers  [12]  and  reported  in  AEDC-TR-71-51  [11]. 

We  shall  proceed  as  follows.  Non-equilibrium  statisti¬ 
cal  mechanics  relates  the  transport  properties  to  time 
correlation  functions  of  the  corresponding  fluxes  [13] . 

Starting  from  these  time  correlation  functions  we  show  in 
Section  II  how  the  first  density  corrections  to  the  transport 
properties  are  determined  by  matrix  elements  involving  a 
triple  collision  operator.  This  triple  collision  operator 
can  be  represented  by  a  binary  collision  expansion.  The 
properties  of  the  binary  collision  operators  are  summarized 
in  Section  III.  In  that  section  we  also  develop  a  diagrammatic 
notation  to  elucidate  the  meaning  of  the  various  terms  in 
the  binary  collision  expansion.  In  Section  IV  we  use  the 
binary  collision  expansion  to  decompose  the  triple  collision 
operator  into  a  series  of  terms  accounting  for  increasingly 
higher  order  dynamical  correlations.  In  Section  V  we  specify 
the  corresponding  three-particle  collision  integrals.  The 
first  density  corrections  n^r  to  the  transport  coeffi¬ 

cients  are  thus  determined  by  a  set  of  collision  integrals 
involving  the  dynamics  of  one,  two,  three  and  four  successive 
collisions  among  three  molecules.  The  results  are  discussed 
in  Section  VI . 
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II .  Transport  coefficients  in  terms  of  collision  operators. 


The  transport  coefficients  for  a  system  of  particles  at 
temperature  T  and  number  density  n  can  be  expressed  as 
integrals  of  time  correlation  functions  of  the  corresponding 
fluxes  [13] 


X  = 


09 

lim  lim  fd t  e"etV"1<J. • J,  (t) > 
3kT  e-*-0  V-*°°  J 


0 

95 


n  =  ihpv  lim  lim  e"etv"1<J  :  J  (t)> 

iUKi  e->o  V->°°  J  n  n 


f 


N 


£•>0  V-*00 


nD  =  i  lim  lim  /  dt  e  e**V  /  v.*v.  (t)  > 

i=l 


(2.1) 


where  k  is  Boltzmann's  constant  and  n  =  <N>/V.  The  average  <  > 
is  taken  over  an  equilibrium  ensemble,  for  which  we  choose 
a  grand  canonical  ensemble  with  temperature  T  =  (k3)  ^  and 
chemical  potential  y.  Thus  for  any  phase  function  X  (x^  ,x^  • .  .x^j) 


=  X (12 . . .  N) 

=  Z  STz  / dxN  TT  4»(vi)W(12...N)X(12...N). 


<x> 


N 


(2.2) 


N 


i=l 


The  grand  partition  function  Z  is  defined  such  that  <1>  =  1; 

2  3/2 

the  activity  z  =  (2irm/3h  )  exp(3y)f  where  h  is  Planck's 

N 

constant  and  m  the  mass  of  a  particle;  dx  =  dx^...dxN  is 
a  volume  element  in  phase  space  and  the  phase  x^  =  (v^,r^) 
represents  the  velocity  v^  and  the  position  r^  of  particle 
i  (i  =  1, 2, . . . ,N) .  The  Maxwell -Boltzmann  factor  ^(v^)  is 
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defined  as 

<Mv)  =  (mB/27T)  3//2exp  (-i-Bmv2)  .  (2.3) 

The  function  W(12...N)  represents  the  statistical  correlations 

between  the  positions  of  the  particles  and  is  given  by 

N 

W(12...N)  =  exp  f-3  )  )  U(r  )  ]  =  ]]  <l+f  )  ,  (2.4) 

j  i<-i 

icj 

where  U(r^)  is  the  pair  potential  as  a  function  of  the 
relative  distance  r. .  =  r.  -  r .  and  f . .  the  Mayer  function. 
For  a  system  of  hard  spheres  with  diameter  o  the  statistical 
correlations  reduce  to  overlap  exclusions  and  the  Mayer 
function  reads 


fij  =  W(ij) -1  =  -0 (d-ri j } , 


(2.5) 


where  0 (x)  =1  for  x  >  0  and  0  (x)  =0  for  x  <  0. 

The  expressions  (2.1)  involve  correlation  functions 
between  the  current  J  at  time  t  and  at  time  t  =  0.  For 

thermal  conductivity  and  viscosity  the  current  consists  of 

K  U 

a  kinetic  term  J  and  a  potential  energy  term  J 


j  =  jk+jd 


(2.6) 


with 


N 


N 


Jx  -  X  W  =  I  <H-§*T)7i  - 

i=i  m 

-  Xw  -  Zm{*iVs  > 


(2.7) 


and 
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N 


N 


5“  - 1  ZZ[u<rij)!-Jijau(i:ij>/aJi3,-Vh0Z?i 

i^j  i=l 


(2.8) 


N 


i^j 


The  symbol  {  }  indicates  a  traceless  symmetrized  tensor,  l 
s 

is  the  unit  tensor  and  h  =  jkT  +  hU  the  enthalpy  per  particle. 
The  self-diffusion  contains  only  a  kinetic  term,  as  shown 
in  (2.1) . 

The  separation  of  the  currents  in  two  parts  (2.6)  implies 
that  X  and  n  can  be  considered  as  the  sum  of  a  number  of 
distinct  contributions 


X  =  xkk+aku+xuk+xuu 


„  -  „KK.KU,UKOlJ 

n  =  n  +n  +n  +n 


(2.9) 


with  obvious  definitions  for  the  individual  terms. 

The  main  text  of  this  report  is  devoted  to  a  study  of 

the  effect  of  triple  collisions  on  the  first  density  correction 

to  the  transport  coefficients.  For  this  purpose  it  is 

KK  KK 

sufficient  to  consider  only  the  kinetic  terms  X  and  n 

The  remaining  terms  are  sometimes  referred  to  as  collisional 

transfer  contributions.  This  terminology  is  slightly  ambiguous, 

since  in  the  original  work  of  Enskog  [7]  collisional  transfer 

*  n  a.  ,  UK  ,UU  .  UK,  UU  .  . 

refers  only  to  X  +X  and  r)  +n  .  The  lowest  order  contri- 

KU  UK  KU  UK 

butions  from  X  +X  and  n  +n  are  proportional  to  the 
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density  n.  They  involve  the  dynamics  of  two  particles 
only  and  are  fully  accounted  for  by  the  theory  of  Enskog. 

In  Appendix  A  we  show  how  these  terms  can  be  derived  from 
a  systematic  density  expansion  of  the  time  correlation 
functions  (2.1).  The  terms  XUU  and  nUtT  account  for  correla¬ 
tions  between  two  pairs  of  molecules;  their  leading  contribu- 

2 

tions  are  proportional  to  n  and  are  not  considered  here. 

In  order  to  discuss  the  time  dependence  of  the  current 

J(t)  in  (2.1)  we  need  to  consider  the  dynamics  of  N  hard 

spheres.  The  dynamics  of  N  interacting  particles  can  be 

tLN 

represented  formally  by  streaming  operators  e  ,  which 

generate  the  solutions  of  the  Hamiltonian  equations  of  motion, 
tL„ 

i.e.  J  (t)  =  e  iNJ(0).  However,  for  hard  spheres  the  streaming 
operators  are  not  defined  when  the  particles  are  overlapping. 
On  the  other  hand  the  streaming  operators  are  always  preceded 
by  the  function  W(12...N)  giving  zero  weight  to  initially 
overlapping  configurations.  One,  therefore,  has  the  freedom 
to  extend  the  definition  domains  of  the  streaming  operators 
to  include  initially  overlapping  configurations.  '  Pseudo 
streaming  operators  e  and  e  ,  defined  for  all  configura¬ 
tions,  which  generate  for  non-overlapping  configurations 
the  actual  trajectories  of  N  interacting  hard  spheres,  have 
been  introduced  by  Ernst,  Dorfman,  Hoegy  and  Van  Leeuwen  [14]  , 
and  will  be  used  in  this  report.  The  operators  LN  and  LN  are 
given  by 
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ln  -  L-+II  Tij  '  Z*  -  L-+ZI 


(2.10) 


where  L0  is  the  Liouville  operator  associated  with  the  free 


streaming  of  the  particles 

N 

L0  =  L0 (12. . .N)  =^vi‘3/9?i 


(2.11) 


and  where  the  operators  and  T\  ^  are  associated  with  a 
binary  collision  between  particles  i  and  j .  The  operators 

and  contain  an  interacting  and  a  non-interacting  term 


T.  .  =  T^.+T?.  ;  T.  .  =  T^.+T*?.  , 

ID  ID  ID  ID  ID  ID 


(2.12) 


such  that 


1  O  f  ^  |  ^  A  i  ■  W  ^ 

Tij  ■  °  id° iD|viD-°ij|S(riJ-aiD>RoJ 


v .  .  •  a .  .  <  0 
ID  ID 


describes  an  interacting  collision, 

T?.  =  -a2  f da.  .  |v. . *a.  . | 6(r .  .-a.  .) 
ID  J  ID  ID  ID  ID  ID 


v .  .  *  a .  .  <  0 
ID  ID 


describes  a  non-interacting  penetrating  collision,  and 

t" .  =  -a2  f da . . j v . . • a .  .  |  6 (r . . -a.  . ) 

ID  'J  ID1  ID  ID1  ID  ID 


v.  .  ‘a.  .  >  0 
ID  ID 


(2.12a) 


(2.12b) 


(2.12c) 


describes  a  non-interacting  separating  collision.  The  symbol 
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A  -  ^ 

0^^  =  cr^j/o  represents  a  unit  vector  (perihelion  vector  of 
the  collision)  and  the'  operator  transforms  the  velocities 


-  13  - 

v^  and  Vj  before  the  collision  into  the  velocities  vt  and  v* 

after  the  collisiont 


-> 

v. 

=  V* 

=  v .  -  (v .  . 

ij 1 

1 

1  13 

V  . 

ij 11 

=  V* 

3 

=  V  .  +  ( V  .  , 

3  13 

t> 

11 

(i^j) 

(2.13) 


where  vij=vi-vj*  Equations  (2.13)  describe  the  "hard  sphere 
interaction"  in  which  molecules  i  and  j  exchange  velocity 
components  along  the  perihelion  vector  . 

To  illustrate  the  meaning  of  these  different  types  of 
binary  collisions  we  consider  in  Fig.  1  the  motion  of  particle 
1  relative  to  particle  2 .  At  a  given  time  t  =  0  particle 
1  is  approaching  particle  2  with  a  relative  velocity  v^2 
and  an  impact  parameter  b12  <  a.  After  a  time  t  =  t^2  the 
center  of  1  will  touch  the  interaction  sphere  of  2  at 
A  (v^2.a12<0)  .  If  the  collision  is  an  interacting  one, 
particle  1  will  proceed  along  its  deflected  path  in  the 
direction  of  v*2.  If  the  collision  is  non-interacting 
penetrating,  particle  1  will  continue  to  move  in  the  direction 
of  v12.  In  the  latter  case  the  particles  will  experience  a 
non-interacting  separating  collision  after  a  time  t  =  t^2, 

•j* 

In  this  report  we  use  forward  streaming  operators  in  contrast 

to  a  previous  publication  on  this  subject  [12] . 
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interaction 


interacting  non-interacting  non-interacting 

collision  penetrating  collision  separating  collision 


Fig.  1.  Geometry  of  a  collision  between  two  hard  spheres. 
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when  particle  1  touches  the  interaction  sphere  of  2  at 
B  (712-S12>0) • 

Since  the  transport  coefficients  (2.1)  are  Laplace 
transforms  of  the  time  correlation  functions,  it  is  convenient 
to  take  the  Laplace  transform  of  the  pseudo-streaming 
operators  (resolvent  operators) 


We  shall  not  indicate  explicitly  that  Le  and  G0  are  N-particle 
operators.  Both  streaming  operators  etLN  and  etLw  or  GN  and 
Gn  describe  the  hard  sphere  dynamics  correctly,  provided 
they  are  used  in  the  appropriate  combination  with  the  overlap 
exclusion.  That  is,  the  operator  etLN  represents  the  dynamics 
correctly  when  written  behind  the  overlap  exclusion  W(12-...N) 
and  the  operator  etLN  when  written  in  front  of  W(12....n). 
Thus,  the  allowed  combinations  are 


W (12 . . .N) e 


W(12.  . .N) 


(2.15a) 


or 

W(12...N)Gn  =  GnW(12...N)  .  (2.15b) 

After  this  review  of  hard  sphere  dynamics,  the  kinetic 
contributions  to  the  transport  coefficients  will  be  related 
to  one-particle  fluxes  J (v)  by  integrating  the  time  correla¬ 
tion  functions  (2.1)  over  the  positions  of  N  particles  and 
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the  velocities  of  N-1  particles  [4] 


KK 


-  -hJ 

3kTz  J 


-t 


dv4>  (v)  (v)  *4^  (v;  0) 


nKK  “  I5kT /a^<v)Jn(?,!%(?!0) 

nD  =  f  dv<|)  (v)  V'Wp  (v;  0) 


(2.16) 


We  represent  the  functions  ¥(v;0)  =  lim  ¥(v;e)  as  the  result 

e->0  ^ 

of  an  operator  T(v;e)  operating  on  the  appropriate  flux  J (v) 

$x(v;e)  =  r(v;e)?x(v)  ,  f  ^  (v;  e)  =  r(v;e)J^(v)  .  (2.17) 

Then  it  follows  from  (2.1)  that 

nHv^e)  =  lim j'dx1*"1  ft  4>  (v^  W(12  . .  .N)  GN  ^  ?  *(2.18) 

V""“  ‘N  "  i=2  j-1 

The  operator  P^  in  (2.18)  is  a  permutation  operator 
which  interchanges  the  labels  of  particles  1  and  j .  The 
corresponding  relation  for  the  self-diffusion  reads 

$D(v;e)  =  r  (v;  e)  v  ,  (2.17c) 

where  T(v;e)  is  obtained  from  (2.18)  by  deleting  the  sum  of 
the  permutation  operators. 

In  order  to  obtain  an  integral  equation  for  the  functions 
¥  in  (2.17)  we  follow  a  method  proposed  by  Zwanzig  [15].  A 
cluster  expansion  of  the  resolvent  operators  leads  to  a  formal 
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density  expansion  for  the  operator  F(v;e).  Next,  we 
determine  the  density  expansion  of  T  \  and  multiply  (2.17) 
with  the  inverse  operator  V  ^  to  obtain  the  desired  integral 
equation. 

For  this  purpose  we  introduce  a  set  of  N-particle  resol¬ 
vent  operators 

n 

GN(12...n)  =  G(12...n)  =  [e-L0  (12  . .  .N) -  T±j  ]  ,  (2.19) 

i  3 

in  which  collisions  (T-operators)  between  only  the  first  n 

particles  (n  <  N)  are  taken  into  account.  Cluster  operators 

are  then  defined  by  the  recursion  relation 

_N 

G  (12 . .  .  n)  )  U  (12 ...  A)  •  (2.20) 

*=1  lCUTc[n] 

The  second  summation  is  taken  over  all  sets  [Z]  of  Z  particles 
which  contain  particle  1  and  are  themselves  contained  in  [n]  = 
(l...n).  For  example, 

G  ( 1)  =  G0  =  {e-L0  (12...N}’"1  =  U(l), 

G (12)  =  U ( 1) +U (12 ) ,  (2.20a) 

G  (123)  =  U(l)+U(12)+U(13)+U(123)  . 

Inversion  of  the  recursion  relation  (2.20)  yields 
Z 

U  (12 . .  .  £)  =  y  ( -1 )£_3C  Y  G  (12  .  .  .  k)  .  (2.21) 

k=l  lC[k]C[Jl] 
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In  particular,  the  first  few  terms  read 
U  (1)  =  G  (1)  =  G0  , 

U  (12)  =  G  (12)  -G0  ,  (2.21a) 

U  (123)  =  G(123)-G(12)-G(13)+G0. 

We  note  that  all  operators  in  (2.20)  and  (2.21)  operate  on 
the  phases  of  N  particles. 

When  we  introduce  this  cluster  expansion  into  expression 
(2.18)  for  T  we  can  start  from  either  W(12...N)G„  or 
GnW(12...N).  In  the  latter  case  we  define  U-operators  by 
replacing  G  with  G  in  (2.20).  Due  to  (2.15b)  both  procedures 
are  equivalent,  but  the  second  representation  is  more 
convenient  in  view  of  the  integrations  over  the  phases  of 
particles  2  and  3 ,  to  be  carried  out  in  the  next  section . 

Next,  the  cluster  decomposition  (2.20)  is  substituted  in 
(2.18)  to  obtain  an  expansion  for  the  operator 

—i  £-1  p.  Z  A, 

r(^l!E)  =  X  71=17  r  I dx*"1  IT  <Mvi)U(12,..«g(12...il)  Y  ?iy  (2.22) 
*-l  1-2  j-I 

where  we  have  introduced  the  equilibrium  distribution  functions 
g (12 .  ..Z)  [16]  : 

N 

dxN_A  [I  4>(vi)WN,  (2.23) 

v^w  N=£  i>=Z+l 

with  g(l)  =  1.  The  reduced  distribution  functions  g(12...A) 
can  be  represented  by  a  power  series  in  the  density.  Here  we 


n  g(12. . .Z)  =  lim 


I  <N-it)!z/ 
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quote  only  those  terms  which  we  shall  need  explicitly  [16]  : 


g  (12)  =  W(12)+nj  dx3<Mv3)W(12)f13f23+*  •*  / 
g  (123)  =  W ( 123 ) + •  •  • 


(2.24) 


The  form  (2.22)  for  the  operator  T  is  appropriate  only  when 
operating  on  functions  J (v)  of  the  velocity  for  which 


J 


dv<f>  (v)  J  (v)  =  0 


(2.25) 


and  this  is  precisely  the  case  for  the  transport  coefficients. 
We  note  that  the  operators  U{12...£)  in  (2.22)  now  operate 
only  on  the  phases  of  £  particles .  The  reason  is  that  the 
terms  operating  on  the  remaining  (N-£)  particles  give  rise 
to  vanishing  surface  terms.  The  operator  r  for  the  self¬ 
diffusion  is  also  given  by  (2.22)  provided  that  we  delete 
again  the  permutation  operators. 

We  expand  the  operator  T  in  a  power  series  in  the 
density 

r  =  r]+nr2+n2r3+- • •  ,  (2.26) 


with 


=  G, 


Pli 


(2.26a) 

(2.26b) 


=  f  dx2<Mv2)U(12)W(12)  £ 

i='l  3 

=  f  dx2dx3<j)(v2)<()(v3)  {|u(123)W(123)+U(12)W(12)f13f23}  ^  P^. 

i=l 

(2.26c) 
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In  order  to  write  in  the  form  (2.26c)  we  have  added  a 
term  U(12)W(12)  ^13^23^13  w^^-ch  vanishes  after  integration 
over  v^  according  to  (2.25).  Again  the  corresponding  terms 
for  the  expansion  of  r  are  obtained  by  deleting  the  permuta¬ 
tion  operators. 

As  pointed  out  by  Zwanzig  [15]  the  density  expansion 

(2.26)  is  not  suitable  for  a  calculation  of  the  functions 

H’tvje)  in  (2.17),  since  the  coefficients  of  the  expansion 

t 

diverge  as  e  in  the  limit  e-*-0.  However,  Zwanzig  has  also 
shown  that  the  inverse  operator  r-1(e)  has  a  finite  limit  as 
e-*0,  at  least  for  the  first  three  terms,  in  the  density  expan¬ 
sion  of  T  ^(e).  This  leads  us,  after  inverting  (2.26)  and 
performing  a  density  expansion,  to  the  following  integral 
equations  for  ¥(v;e): 

J(v)  =  limr-1(v;e) ¥(v;e)  =  -{nl, (v) +n2I^ (v) +...}¥ (v; 0) ,  (2.27) 

e->0  Zi 

with 

I  (v)  =  lirnGo^T- (v; ejGo1,  (2.28a) 

e->0 

I o ( v)  =  lim{Gr1r. (v;£)G;1-G;1r9(v;E)Go1r9(v;£)G;1>.  (2.28b) 

J  E+0  J  Z  Z 

The  operator  G0 ^  can  be  replaced  with  e  when  acting  on 

functions  of  the  velocities  alone.  In  deriving  (2.27)  we 

have  used  the  property  that  limG7"^  (v; e)  =  lime'i' (v;  e)  =  0. 

^  e+0  e+0 

The  operator  ^(v)  i-s  the  Boltzmann  collision  operator 
which  will  be  reduced  to  its  familiar  form  in  the  next 
section;  the  operator  I^(v)  is  a  triple  collision  operator. 

In  the  integral  equations  (2.27)  both  operators  act  on  a 
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function  of  the  velocity  alone.  For  the  self-diffusion 
we  obtain  similarly  the  integral  equation 

v  =  -Cnl2(v)+n2i3(v)  +  -  •  -}^D(v;0)  ,  (2.29) 

where  the  operators  Ik(v)  are  given  by  (2.28)  after  replacing 
rk  with  rk. 

The  solution  ¥(v;0)  of  the  integral  equation  (2.27)  can 
be  written  in  the  form  of  a  density  expansion 

¥ (v; 0)  =  -n"1{l"1-nl“1I3l21+***}J(v)  ,  (2.30) 


and  a  similar  result  follows  from  (2.29).  When  the  formal 
solution  (2.30)  is  inserted  into  (2.16),  we  obtain  a  density 
expansion  for  the  kinetic  parts  of  the  transport  coefficients 


,KK  ,  , KK  „ 

A  =  A0+A^  n+* • 


KK  _  KK 

n  =  ho+n-L  n+‘ 


(2.31) 


nD  =  Dq+D^n+*** 

The  lowest  density  contributions  are  in  general  of  the  form 


-/< 


-1 


dv<}>  (v)  J  (v)  1 2  (v)J(v) 


(2.32a) 


while  the  first  order  density  corrections  have  the  form 

/dv<J>  (v)  J  (v)  l"1I3l21J  (v) ,  (2 . 32b) 


as  follows  from  (2.16)  and  (2.30). 
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In  order  to  express  the  quantities  (2.32)  as  matrix 

elements  of  collision  operators  we  introduce  the  functions 
4  -± 

A,  B  and  C,  which  are  determined  by  the  integral  equations 
I2(v)A(v)  =  -Jx  (v)  # 

I2(v)B(v)  =  -Jn(v)  ,  (2.33) 

I2(v)5(v)  =  -v 


From  (2.32) ,  (2.33)  and  the  symmetry  of  the  Boltzmann  collision 
operators  I2,  it  follows  then  that  the  low  density  values 
A0/  no  and  D0  are  given  by 

Ac  =  ~  o"  lA'Al  (2)  * 

3kT: 

Ho  =  [B#B]  (2)  ,  (2.34) 

D0  =  “j  t£,C]  (2)  . 

KK  KK 

For  the  first  density  corrections  A^  ,  n-^  and  we  obtain 

A«  =  1  [ifft  O), 

1  3kT^ 

KK  _  JL _  (3) 

^1  lOkT  IB/BJ  •  (2.35) 

D1  =  J  [C,6]  • 

For  convenience  we  have  introduced  a  short  hand  notation 
defining  matrix  elements  associated  with  the  operators  1^  and 
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[f,g]  (k)  =  J dv<ji  (v)f  (v)*Xk(v)g(v) 
[f  ,g]  dv<ji  (v)  f  (v)  *Ik(v)g  (v) 


(2.36a) 

(2.36b) 


where  g(v)  and  f (v)  are  vector  or  tensor  functions  of  the 
velocity  v  and  *  indicates  the  appropriate  scalar  product. 


21 


AEDC -TR-72-142 


III .  Dynamics  of  hard  spheres 

Once  the  coefficients  in  the  expansion  (1.1)  for  the 
transport  properties  are  related  to  collision  operators, 
the  problem  is  reduced  to  an  analysis  of  the  dynamics  of 
two,  three,  etc.  isolated  hard  spheres.  To  describe  the 
dynamics  we  introduced  binary  collision  operators  in  (2.12). 
We  now  summarize  some  of  the  properties  of  the  T-operators 
and  elucidate  how  they  enable  us  to  give  a  compact  and 
systematic  description  of  collision  sequences. 

The  T-operators  correspond  to  binary  collisions  which 
may  be  interacting  or  non-interacting.  The  G0 -operators , 
defined  in  (2.14),  correspond  to  free  streaming  of  the 
particles.  Products  of  T-operators  and  G0 -operators  corre¬ 
spond  to  sequences  of  successive  binary  collisions  among  the 
particles,  and  can  be  represented  by  diagrams.  The  G0- 
operators  are  represented  by  labeled  particle  lines  (labels 
1,2,3)  indicating  the  free  trajectories  of  the  particles 
between  collisions.  We  adopt  the  convention  that  time 
increases  when  the  diagrams  are  read  from  bottom  to  top  as 
indicated  by  arrows  on  the  particle  lines.  The  T-operators 
correspond  to  the  basic  binary  collision  events  and  are 
represented  schematically  in  Fig.  2.  The  collisions  are 
assumed  to  occur  at  a  given  time  x  =  x0  and  in  Fig  2a  we 
indicate  that  the  velocities  change  in  an  interacting 
collision.  In  Fig.  2b  and  2c  the  velocities  do  not  change. 
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Fig.  2.  Schematic  representation  of  (a)  an  interacting  collision,  (b)  a  non¬ 
interacting  penetrating  collision,  (c)  a  non-interacting  separating 
collision,  and  (d)  the  condition  that  the  two  particles  overlap  at 
time  t=Tq. 
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but  we  place  marks  just  above/below  the  collision  time  in  the 
non-interacting  penetrating/separating  collision  to  indicate 
that  the  particles  overlap  just  after/before  the  collision. 
Occasionally,  it  will  be  necessary  to  require  that  the  parti¬ 
cles  overlap  at  a  given  time  x  =  x0;  we  do  this  by  connecting 
their  trajectories  with  marks  as  shown  in  Fig.  2d. 

Using  the  foregoing  symbolic  representation,  we  can 
construct  diagrams  representing  sequences  of  collisions. 

As  examples  we  show  in  Fig.  3  the  diagrams  associated  with 
the  various  terms  in  the  operator  products  Ti2G°T13  and 
T12G0T13.  The  operator  T12  is  only  different  from  zero 
when  particles  1  and  2  are  in  contact.  The  operator  T^3 
requires  that  at  the  same  time  the  phase  of  particle  3  is 
such  that  the  collision  between  1  and  2  is  followed  by  a 
collision  between  1  and  3,  i.e.  x13>x12>  For  interactin5 
and  non-interacting  penetrating  collisions  the  time  ordering 
refers  to  xa  and  for  non-interacting  separating  collisions  to 
x+,  where  a  indicates  the  appropriate  pair.  The  operators 
act  on  functions  of  the  positions  and  velocities  of  the 
three  molecules.  The  effect  of  the  operator  product  is  to 
transform  the  positions  and  velocities  of  the  molecules  just 
prior  to  the  first  collision  at  the  bottom  of  the  diagram 
into  those  just  after  the  last  collision  at  the  top  of  the 
diagram^.  As  mentioned  earlier,  any  penetrating 

collision  will  be  followed  by  a  separating  collision  at  a 


^Notice  that  the  bottom/top  of  the  diagram  is  determined  by  the 
first/last  collision.  The  extensions  of  the  particle  lines  below 
the  first  collision  and  above  the  last  collision  are  shown  for 
convenience  only;  they  do  not  represent  G0 -operators . 
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later  time.  However,  the  diagrams  are  non-committal  about 
the  number  of  and  the  time  at  which  additional  non-inter- 
acting  collisions,  not  indicated  explicitly,  will  occur. 

For  example?  in  Figs.  3c  and  3d  a  separating  collision  between 
1  and  2  will  occur  at  x  =  x^ •  but  x^  ma^  e^ther  smaller 
or  larger  that  x^*  As  mother  example,  we  show  in  Fig.  4 
the  diagrams  associated  with  the  two  terms  of  the  operator 
product  f23T12G°T13‘  In  these  collision  sequences  particles 
1  and  2  are  colliding,  while  2  overlaps  with  3.  These 
diagrams  can  roughly  be  thought  of  as  space-time  plots  of 
the  centers  of  the  molecules,  with  the  time-axis  vertical  and 
the  (one-dimensional)  space-axis  horizontal,  and  they  allow 
a  direct  geometrical  interpretation  of  the  collision  operator 
products . 

The  operators  satisfy  a  number  of  properties  listed  in 
Table  I.  Proofs  of  these  relations  are  either  given  by 
Ernst  et  al.  [14]  or  are  directly  implied  by  their  results. 

The  relations  (1.1)  through  (1.6)  have  an  obvious  dynamical 
meaning,  when  interpreted  geometrically.  In  (1.8)  the 
operation  of  time  reversal  is  represented  by  the  operator  T 
which  reverses  the  directions  of  the  velocities  of  all 
particles,  while  their  positions  remain  unchanged. 

Any  product  of  T-  and  G0 -operators  can  be  interpreted 
as  a  specific  sequence  of  collisions.  However,  some  sequences 
may  not  be  allowed  according  to  the  laws  of  mechanics . 
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TABLE  I 

PROPERTIES  OF  BINARY  COLLISION  OPERATORS 


faTa  =  0  '  V.  =  0  '  faG-Ta  =  0  '  TaG»fa  =  0 

TaG«Ta  =  0  '  *ctG»Tct  =  0  ,  TaG0Ta  =  0 

These  relations  apply  also  to  i-  and  n-collisions 

separately.  In  these  relations  G0  may  be  replaced  by 

G0Tn  G0 . .  .T**  G0  if  a.  a  for  i  =  l,2...Jt;in  addition  any 
al  aZ  1  _ 

(a^a)  may  be  replaced  by 


faTB  =  Va'  faT(J  "  V« 


faG# 

+ 

G«T„G»  = 

G0fa  +  G8TaG0 

faG° 

+ 

GoT^C.  = 

G0fQ  +  G0t£g0 

faG° 

= 

-faG"fa 

-  faGoTjG0 

G°fa 

- 

~foG° fa 

-  GoT^.fa 

TjCo 

- 

-TnG-f 
ot  a 

-  t;g0t^Go 

G0T^ 

= 

-f  G0f" 
a  a 

_  G0TnG0Tn 

P  1  ■>' 

v .  +  v . 


5  2  2 

vf  +  VT 

L1  1 


f  1  - 

V.  +  V  . 


.  T  V  .  p  v .  +  V  .  _ 

1  3  =  0,  /dv.dv ,  1  3  T  f(x. ,x  )  =0 

2  2  J  1  3  2  2  1J  1  J 


v.  +  v . 
1 


fi* ijTijf(*i'fy  -  J^ijTijf<Vi,Vj> 
where  f(x^,Xj)  and  f(vi#y..)  are  arbitrary  functions. 


=  Tt  T-1  =  Tt  T  _1 

a  a  a  a 

The  dagger  represents  hermitian  conjugation;!  is  the  opera¬ 
tion  of  time  reversal. 
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Of  course,  the  laws  of  mechanics  are  here  extended  to  include 
the  possibility  of  non-interacting  collisions.  The  defini¬ 
tion  of  the  T-operators  implies  that  any  operator  product 
associated  with  a  forbidden  collision  sequence  vanishes 
automatically.  Some  of  these  forbidden  collision  sequences 
are  of  a  trivial  nature  and  follow  from  the  properties  of 
Ta  and  Tq  compiled  in  Table  I.  For  example,  it  follows  from 
(1.1)  that 

TxG«T6G“Tc1  -  0  .  (3.1) 

This  relation  expresses  the  fact  that  a  pair  of  particles 
cannot  recollide  after  a  collision,  unless  the  trajectory  of 
at  least  one  of  the  two  particles  is  deflected  by  an  interme¬ 
diate  interacting  collision  with  the  third  particle. 

In  addition  to  such  obvious  restrictions,  there  are 
several  restrictions  which  are  not  of  a  trivial  nature,  but 
are  consequences  of  theorems  concerning  the  dynamics  of 
three  hard  spheres.  In  Table  II,  we  present  lemmas  that 
express  forbidden  collision  sequences.  These  lemmas  are 
consequences  of  the  following  two  theorems: 

THEOREM  1 .  Three  hard  spheres  (with  equal  masses  and 
diameters)  cannot  undergo  more  than  four  successive  inter¬ 
acting  collisions.  Moreover,  any  sequence  of  four  successive 
interacting  collisions  is  either  of  the  type  described  by 
the  diagram  of  Fig.  5,  or  else  obtainable  therefrom  by  inter- 
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TABLE  II 


LEMMAS  INDICATING  FORBIDDEN  COLLISION  SEQUENCES 


The  lemmas  hold  for  the  interacting  and  non-interacting  terms 
of  each  T  -  or  T  -operator  separately.  The  lemmas  remain 
valid  upon  insertion  of  any  number  of  additional  Tn  -  or  T11  - 
operators.  Moreover,  T^  may  be  replaced  by  T^. 
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changing  particle  lables  and/or  reversing  the  direction  of 
time. 

THEOREM  2 .  In  the  recollision  sequence  defined  by  the 
diagram  of  Fig.  6,  it  is  not  possible  for  spheres  2  and  3 
to  collide  or  overlap  at  any  intermediate  time:  i.e. 
r32(f)>a  for  Tr5-T-Tiii  • 

Theorem  1  was  stated  by  Sandri  and  coworkers  [17]  and 
proved  by  Murphy  and  Cohen  [18]  .  Theorem  2  was  proved  by 
Hoegy  and  Sengers  [12].  For  further  details  the  reader  is 
referred  to  AEDC-TR-71-51  [11]. 

In  addition  to  the  concepts  already  used  in  describing 
the  mechanics  of  hard  sphere  collisions,  it  is  convenient 
to  introduce  the  notion  of  "overlap  collisions",  in  order 
to  distinguish  between  statistical  and  dynamical  correlations. 
An  overlap  collision  is  a  collision  of  any  type  (i.e.  inter¬ 
acting,  penetrating  or  separating)  between  two  particles 
which  occurs  while  at  least  one  of  the  two  colliding  particles 
overlaps  with  the  third  particle.  The  collision  between  1 
and  2  in  Fig.  4  is  an  example  of  such  an  overlap  collision. 

The  collision  between  1  and  3  in  Figs.  3c  and  3d  would  be 
an  overlap  collision  in  the  case  Ti2>Ti3*  ^e  may  father 
distinguish  between  single-overlap  collisions  and  double¬ 
overlap  collisions.  In  a  single-overlap  collision  only  one 
of  the  two  colliding  particles  overlaps  with  the  third 
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Fig.  6.  Recollision  sequence  for  Theorem  2. 
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particle  and  in  a  double  overlap  collision  both  colliding 
particles  overlap  with  the  third  particle. 
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IV..  Expansion  of  the  collision  operators. 


The  coefficients  of  the  density  expansion  for  the 
transport  properties  are  determined  by  matrix  elements  or 
collision  integrals  (2.34)  and  (2.35)  which  involve  the 
collision  operators  and  given  by  (2.28).  In  this 
section  we  show  how  these  collision  operators  can  be  de¬ 
composed  into  terms  that  correspond  to  collision  sequences 
involving  increasing  numbers  of  successive  binary  collisions. 

For  this  purpose  we  recall  that  the  collision  operators 
I2  and  I^  are  related  to  the  resolvent  operators  G(12) 
and  G(123)  via  (2.26)  and  (2.21).  From  (2.19)  it  follows 
that  the  resolvent  operators  satisfy  the  integral  equations 


G(12. . .A) 


G„+G, 


l 


T  G (12 . 
a 


•  A) 


f 


a 

where  a  runs  over  all  pairs  among  the  Z  particles. 

collision  expansion  is  obtained  upon  interation  of 
00 


G (12 . . . Z) 


s-0  a 
In  particular 


(4.1) 


A  binary 
(4.1) 


(4.2) 


G  (12 )  =  G0+G0T12G0  .  (4.3) 

In  equations  (4.1)  -  (4.3)  the  operators  G  and  T  may  be  replaced 
by  G  and  T,  respectively.  We  also  note  that  those  terms  in 
(4.2)  with  two  successive  T-operators  with  the  same  index  a 
vanish  automatically  as  a  result  of  (1.1).  This  is  also 
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the  reason  why  (4.3}  terminates  after  the  second  term. 

We  first  consider  the  Boltzmann  collision  operator 
I2 (v) .  If  we  now  substitute  equation  (4.3)  for  G(12) 
into  the  expressions  (2.28a)  for  I2(v),  using  (2.26b)  and 
(2.21a),  we  obtain 


=  f  dx2<j>(v2)T12(l+P12)  = 

/  dv2*(v2)o  J  da12 I V12 ' ^12  I  (Ra19“1)  (1+P12)  ' 


(4.4) 


V12*°12  <  0 


which  is  the  familiar  form  of  the  Boltzmann  collision  opera¬ 
tor.  The  operator  I2(v3)  for  the  self-diffusion  is  obtained 
from  (4.4)  by  deleting  the  permutation  operator  P12* 

In  order  to  derive  a  similar  expression  for  the  triple 
collision  operator  I3  (v)  it  is  convenient  to  symmetrize  over 
the  labels  2  and  3  in  (2.26c) .  If  we  then  inseft  (2.21)  and 
iterate  (4.1)  once,  we  obtain  from  (2.28b)  after  rearranging 
some  terms 


I3(vl> 


with 

T  (123) 


j 

=  |  jTdx2dx3<()  ( v 2 )  <j>  (v3)T  (123)  ^  V1L 


(4.5) 


i=l 


=  lim  )t  G(123)W(123)- }  T  G0  + 
e+0  al  ^  ai 

al  al 

-H  v-S^V-V-’K 


(4.6) 


°1^°2 
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The  indices  and  a2  in  the  summations  run  over  the  pairs 

(12) ,  (13)  and  (23).  To  obtain  a  symmetric  expression  we 

have  added  terms  starting  with  T23  and  terms  ending  with 

T. .  p  (k^i,j);  they  vanish  as  a  result  of  (1.7).  For  the 
i]  lk  J 

self-diffusion  one  derives  similarly 


h  (vi> 


dx2dx34)  (v2)<j>  (v3)T(123)  f 


(4.7) 


where  T(123)  is  again  given  by  (4.6).  As  a  next  step  we 
replace  G(123)W(123)  by  W(123)G(123)  according  to  (2.15b) 
for  N=3  and  insert  the  expansion  (4.2)  for  G(123).  After 
rearranging  terms  and  using  (1.1)  and  (1.4)  we  find  [12] 


T (123) 


a. 


“l*»2 


Va,G«V 


(( x)  (a1jiy^a2 ) 


(4.8) 


+  ZI  a+f  )T  ,G.T  Y[Z  G.*a. 

°1*»2  1  2  i-1  °i 


(a1^a2) 

Some  of  the  conditions  on  the  summations  in  (4.8)  are  auto¬ 
matically  satisfied  due  to  (1.1).  It  is  always  understood 
that  the  limit  e-»-0  has  to  be  taken  at  the  end  of  the  calcula¬ 
tions,  but  we  omit  indicating  this  limit  explicitly. 

Expansion  (4.8)  still  contains  many  terms  corresponding 
to  collision  sequences  that  are  not  allowed  by  the  laws  of 
mechanics.  In  the  last  term  of  the  right  hand  side  of  (4.8) 
the  term  with  f  vanishes  as  a  result  of  lemmas  (II. 1)  and 
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(II. 2).  This  means  that  the  conditions  imposed  by  the 
presence  of  (1+fy)  are  automatically  satisfied.  Further¬ 
more,  lemmas  (II.4)  -  (II. 7)  imply  that  the  expansion 
terminates  after  products  of  four  T-operators.  Thus 

TU23)  =££B£Y 
a 

*  Y1  (4'9) 

a  3 

+  [TaG.T£G.T0G„TY+faGoTeGsT£GoT6]  , 

a  3 

where,  from  now  on,  we  use  the  convention  a^3^Y^a*  In  (4.9) 
we  have  also  omitted  those  combinations  of  three  or  four 
T-operators  which  vanish  according  to  (1.1). 

As  mentioned  in  the  introduction  we  want  to  make  a 
distinction  between  statistical  correlations  and  dynamical 
correlations.  Statistical  correlations  are  determined  by 
those  points  in  configuration  space  for  which  two  or  more 
particles  are  overlapping.  In  the  expansion  of  the  collision 
operators  they  lead  to  the  presence  of  "overlap  collisions". 
Therefore,  we  want  to  separate  the  collision  sequences  into 
collision  sequences  with  and  without  overlap  collisions. 

The  expansion  (4.9)  ,  however,  is  not  yet  in  a  form  which 
makes  this  distinction  evident.  The  reason  is  that  only 
those  configurations  for  which  the  particles  overlap  at 
the  time  of  the  first  collision  are  indicated  explicitly  in 
(4.9)  by  Mayer  functions.  In  addition,  however,  the  collision 


T  + 
a 


I? 


W-V 
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sequences  may  also  generate  overlap  collisions  at  a  later 
time,  not  indicated  explicitly  in  (4.9).  Such  overlap 
collisions  will  occur,  if  a  particle  has  penetrated  a 
second  particle  and  collides  with  a  third  particle  before 
it  separates  from  the  second  particle.  Thus  we  consider  in 
more  detail  those  collision  sequences  in  (4.9)  which  contain 
an  intermediate  Tn-  operator.  The  formal  analysis  is  most 
easily  carried  out  by  using  (1.6),  which  leads  to  the  follow¬ 
ing  decompositions. 


TaG-TSG"TV  - 


„n 


n„  =n. 


TaGoT;G0T“G0TY  «  ^G.TJG.T^.TJG.T  ^0G.T;GlT;Gtf0T. 


(4.10) 

(4.11a) 


TaG0T^G0T;GoTe  = 


-t..g0t^g0t^g0t1g0t0-t  G-T^G 


a 


3  °  3 


a 


rir 


T1G„T, 


(4.11b) 


The  tertns  in  (4.11)  represent  sequences  of  four  successive 
collisions.  The  last  term  in  (4.11a)  vanishes  due  to  lemma 
(ll.3)  and  the  last  term  in  (4.11b)  vanishes  as  a  result  of 
(1.1).  Since  the  remaining  terms  on  the  right  hand  side 
of  (4.11)  do  not  contain  any  Mayer  functions,  we  thus  conclude 
that  sequences  of  four  successive  collisions  cannot  lead  to 
any  overlap  collisions  [11,19].  If  we  substitute  (4.10)  and 
(4.11)  into  (4.9)  we  obtain  the  following  expansion  of  the 
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collision  operator 


T (123)  =  Y  T  (123} 


P=1 


with 


(123)  =  £  Tg(o)  =  Y  £6fTT«  ' 
a  a 

T2(123)  =  yy  Tg  (a  /  3)  i  yy  VfYGo“GaT?G*VTe 
a  $  a  3 

T3(123)  =  yy  [TR(a,e)+Tc<a,B)+TH(a,|3)  J  = 


a  p 


T, 


<123>  =  S 


Trc  (a ,  p)  +Tcr  (a  ,  p)  +Trh  (a ,  p)  +TRR  (a ,  p)  1 


E  IE  [TaG„TiGoTjG„VT(iG„TjGoT£G„Te 

a  p 


-TaGoT13Go^G0^G0TY-TaG0T2G0^G0T^GoT3  j  . 


(4.12) 

(4.13) 

(4.14) 

(4.15) 


(4.16) 


The  first  term  T_(a)  =  fQf  T  requires  that  pair  a  collides 

r»  p  y  a 

while  both  molecules  are  overlapping  with  the  third  molecule. 
It  contains  the  dynamics  of  only  one  collision,  which  is 
evidently  a  "double  overlap"  collision.  This  term  is  pre¬ 
cisely  the  term  given  by  the  theory  of  Enskog  as  will  be 
demonstrated  in  the  next  section.  We  therefore  refer  to 
this  term  as  the  double  overlap  or  Enskog  term. 

To  elucidate  the  meaning  of  the  second  term  it  is 
advantageous  to  write  it  as  the  sum  of  three  terms,  using 
(  1.5)  , 
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Tg  (a ,  3)  =  Tss<a,3)+TSN(a,3)+TNS(a,3)  = 

=  -Ta{fYG0fY+fYG0T^G0+G0T^G0fY}Te  (4.17) 

This  term  contains  "the  dynamics  of  two  successive  collisions 
involving  respectively  pair  a  and  pair  3.  In  the  term 
Tss  (a,  3)  the  first  and  last  collisions  are  single  overlap 
collisions,  while  TgN(a,3)  or  TNg(a,3)  corresponds  to  the 
case  where  the  first  collision  is  a  single  overlap  collision 
and  the  last  collision  is  a  non-overlapping  collision,  or 
vice  versa.  We  therefore  refer  to  term  T  (a,3)  as  the 
single  overlap  term. 

The  terms  in  Tg(123),  given  by  (4.15),  represent 
sequences  of  three  successive  collisions  to  which  we  refer 
as  recollisions  ,  cyclic  collisions  and  hypothetical  collisions. 

The  terms  in  T4(123)  represent  sequences  of  four  successive 
collisions,  which  are  combinations  of  recollisions,  cyclic 
collisions  and  hypothetical  collisions.  We  emphasize-  that 
both  T3(123)  and  T4(123)  do  not  contain  any  overlap  collisions. 
The  geometrical  meaning  of  the  various  collision  sequences 
will  be  elucidated  in  the  next  section. 

The  triple  collision  operator  (v) ,  which  is  related 
to  T(123)  via  (4.5),  is  a  symmetric  operator  [11,12).  To 
exhibit  the  symmetry  we  take  the  hermitian  conjugate  of  the 
operator  T(123).  Using  (1.8)  and  the  relation  G^  =  Tg0T  ^ 
we  find  for  the  Enskog  term 
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=  fflf  TT  T-1  , 

E  6  y  a 

and  after  appropriate  relabeling  of  the  particles 
TgS(a,B)  =  TTss(a,f3)T"1 
T+N(a,B)  =  TT^a^JT"1 

Tj(a,B)  =  TTR(a,3)T_1 
T<t(a,B)  =  TTC  (a ,  6)  T-1 

Tj(a,B)  =  TTjjCa^)!-1 

T^c(a,3)  =  rTCR(a,B)T-1 
TRH(a,B)  =  TTHR(a,B)T-1 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


Thus  the  triple  collision  operator  I3 (v)  is  symmetric  in  the 
space  of  functions  f (v)  of  the  velocity  which  have  a  definite 
parity  (±)  under  the  operation  of  time  reversal,  i.e. 

Tf(v)  =  f (-v)  =  ±f(v).  (4.22) 


The  functions  of  interest  in  calculating  matrix  elements 
of  above  operators  (4 .18)  -  (4 . 21)  are  the  solutions  of  the  linear 
inhomogeneous  integral  equations  (2.33).  These  solutions 
have  indeed  a  definite  parity,  while  the  linearized  Boltzmann 
operator  ^(v)  is  isotropic  in  velocity  space.  The  isotropy 
of  I3(v)  in  velocity  space  ensures  that  (v)  does  not  have 


42 


AEDC-TR-72-142 


any  matrix  elements  between  functions  of  different  parity. 
A  proof  of  the  symmetry  of  the  Zwanzig  form  of  the  triple 
collision  operator  (v)  for  a  general  potential  has  been 
given  by  Ernst  [20  ]  . 
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V.  Specification  of  collision  integrals 


The  expansion  (4.12)  for  T(123)  implies  an  expansion 
for  the  first  density  correction  (2.35)  to  the  transport 


coefficients 


X 


>?-  E 


’ly  lOkT 


I  ^ 


(5.1) 


Dx  ■  IDiu  -  i  I  ^ 


y=l  y=l 

where  the  matrix  elements  are  obtained  from  (2.36)  and  (4.5)  , 
if  T(i23)  is  replaced  by  T^(123)  as  given  in  (4.12)  -  (4.16). 
The  functions  A(v),  B(v)  and  C  (v)  ,  determined  by  the  Boltzmann 
equation  (2.33),  are  usually  given  in  the  form  of  a  Sonine 
polynomial  expansion  [1] .  A  numerical  representation  of  these 
functions  has  been  obtained  by  Brooker  and  Green  [21]  .  In 
general,  a  calculation  of  the  transport  coefficients  (5.1) 
requires  the  evaluation  of  a  set  of  triple  collision  matrix 
elements  of  the  form 

WfXljJ3*  ■  j  y'dv1dx2dx3$(123)^(vJL)*Ty  (123)  ^  x(v±)  , 


[4»rXl^3>  =  jf  dv1dx2dx3$(123)^(v1)*Ty(123)X(v1) 


(5.2) 


where  we  have  introduced  the  shorthand  notation 


$(123)  =  4>  (vx)  <|)  (v2)  4>  (v3)  . 


(5.3) 
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It  is  convenient  to  symmetrize  these  integrals  by  inter¬ 
changing  the  velocity  variables  to  yield 

3  3 


EU> 


'Xly3)  =  |/ dv1dx2dx3$(12  3)  £  (v^  *T^  {123 )  £  X(Vj)  , 


i=l 


j  =  l 


(5.4) 


|  f  dv1dx2dx3$  (123)  ^  *Ty  d23)  X  (v±) 

i=l 

In  this  section  we  shall  investigate  the  structure  of 
these  collision  integrals  and  elucidate  how  they  are  related 
to  specific  collision  sequences.  The  collision  integrals 
corresponding  to  the  first  term  in  the  expansion  (5.1)  are 
obtained  by  substituting  (4.13)  into  (5.4).  Since  the 
three  terms  of  (4.13)  yield  equal  contributions,  we  may  write 

2  2 

■  I  /dVx2dX3f13£23$<123>  X,,(^i)*fl2Z'X(V  •  <5-5) 

i-1  j-1 

where  the  terms  i=3  and  j  =  3  vanish  as  a  result  of  (1.7)  .  We 
note  that  the  operator  T^2  contains  a  5-function  which  restricts 
the  distance  r^2  to  the  value  o.  Using  the  relation 


[/■ 


drl3f13f23 


r12=o 


5  3 

12  710  ' 


(5.6) 


which  is  the  value  of  the  first  density  correction  to  the 
radial  distribution  function  g(r12)  at  ri2  =  we  conclude 

[ipfXli3)  =  if  7103 I^fXl  r 


(5.7) 


(2) 

where  [ ^ , x ]  is  the  well  known  binary  collision  integral 
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associated  with  the  Boltzmann  equation.  For  the  self-dif fu¬ 
sion  we  find  similarly 


[*,X] 


(3)  = 
1 


_5 

12 


na3  [t|»#X]  (2) 


(5.8) 


It  thus  follows  from  (5.1)  and  (2.34)  that 

.  5  3. 

X11  "  "12  710  X° 

5  3 

"Y2  9 

Dn  =  “if  7ra3°° 


(5.9) 


which  are  precisely  the  contributions  predicted  by  the  theory 
of  Enskog  [7] . 

We  now  turn  our  attention  to  the  contributions  for 
y= 2,3,4,  not  contained  in  Enskog' s  theory.  Since  the  collision 
integrals  (5.4)  are  symmetric  with  respect  to  permutations  of 
the  particle  labels,  it  suffices  to  evaluate  operators  T^(123) 
for  one  permutation  only.  Thus 


lY-Xl,!3’  = 

J1  dv1dx2dx3$ (123) 

C*,xl‘3>  - 

J'  dv^dx2dx3$  (123) 

f  *<?1).TJ  (12,13)  Y  X(v.) 
ijl  J-l 


^  ip  (vi)  *T^  (12,13)  x(vi) 


i=l 


where 


(5.10) 


T2 (12,13)  =  Tss(12,13)+Tsn(12,13)+Tns (12,13) 

T3(12,13)  =  TR(12,13)+TC(12,13)+TH(12,13)  ,  (5.11) 

T4 (12 ,13)  =  Trc(12,13)+Tcr(12,13)+Trh(12,13)+Thr(12,13)  , 
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as  defined  by  (4 .14) - (4 . 17) .  Note  that  in  (5.10)  the  term 
i=3  vanishes  as  a  result  of  (1.7).  It  should  be  noted 
also  that  the  summation  over  the  particles  j  in  (5.10) 
contains  only  those  two  particles  which  are  involved  in 
the  last  collision. 

The  matrix  elements  in  (5.10)  contain  binary  collision 
operators  (2.12)  and  describe  sequences  of  binary  collisions. 

In  this  section  the  collision  sequences  corresponding  to  each 
matrix  element  will  be  represented  by  diagrams,  using  the 
notation  developed  in  section  III,  the  only  difference  being 
that  we  don't  find  it  necessary  to  specify  whether  the  first 
and  last  collisions  are  interacting  or  non-interacting.  In 
addition  the  expressions  for  the  matrix  elements,  which 
still  contain  6-functions  through  (2.12),  are  reduced  to 
fourteen  dimensional  integrals  over  regions  of  phase  space, 
where  the  corresponding  collision  sequence  is  geometrically 
possible.  Furthermore,  we  specify  the  integrands  completely 
in  terms  of  the  integrating  variables.  Since  the  reduction 
of  the  matrix  elements  is  rather  involved,  we  give  the 
details  for  two  specific  examples  in  a  separate  appendix  B, 
and  simply  list  the  results  here. 

We  consider  first  the  single  overlap  term  u=2,  which 
can  be  written  as  a  sum  of  three  terms  according  to  (4.17) 
and  (5.10) 

[*,X]j3)  =  ^rX]^s)  +  [^'x]SN+[,f''x]NS)  '  (5.12) 
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where 


[*»Xlgs>  "  -/  dv1dx2dX3$(123)^(v1/v2)*T12f23G0f23T13x(v1,v3)  , 


t^rXl 


t^fXl 


(3)  =  _ 


SN 


(3)  =  _ 


NS 


(5.13a) 

/  aviax2dX3t(123)*(v1,v2)»f12f23G.^3G0T13x(v1,V3), 

(5.13b) 

J  dv1dx2dx3$(123)ii/  (^^^,^2)  23T13x (v1#v3)  , 

(5.13c) 


and  where  we  have  introduced  the  shorthand  notation 


*(Vi,Vj)  =  ip  (v±)+^  (Vj  )  , 


(5.14) 


X  (v±  # vj )  =  x(vi)+x(Vj)  . 


The  terms  in  (5.13)  are  represented  by  the  diagrams  in  Fig.  7. 
In  order  to  express  the  matrix  elements  in  a  compact  form 
we  define 


Ao . .^(^i'^j)  =  + <vi* *vj*) -V ^vi'vj>  f 

where  the  velocities  v. *  and  v.*  are  functions  of  o- v. 

1  J  ^-]  i 

and  Vj  as  given  in  (2.13).  Using  the  results  of  appendix  B 
we  find  for  the  SS-collisicn  sequence 


(5.15) 


(3) 

ss 


=  -  fan*  (123)  Aa  ip(v1,v2)*Aa  x(v1,v3)  ,  (5.16) 

12  A  ‘  13  i  J 


(SS) 


where  the  volume  element  is 
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dQ  =  02 | v12‘012 |d012dr13dv1dv2dv3  .  (5.17) 

The  conditions  determining  the  (SS) -integration  region 
are  given  by  the  diagram  in  Fig.  7a,  and  read  explicitly 

V12-S12  >0  ,  r23  <  a  '  <5‘18a) 

(r13xv13)2-(ov13)2  <  0  :  t3  >  0  ;  (5.18b) 

t1  <  t2  .  (5.18c) 


The  integrating  variables  represent  the  positions  and  veloci¬ 
ties  just  after  the  first  collision.  That  is,  r12  =  o^2 
and  r^3  refer  to  the  relative  positions  at  the  time  of  the 

first  collision  at  the  bottom  of  the  diagram  and  (v)  = 

B 

vl'v2'v3  to  the  velocities  between  the  first  and  the  last 
collision. 

The  other  symbols  are  defined  in  terms  of  the  integrating 
variables  as  follows 


13  = 

r13+v13Tl 

• 

t 

r23 

4-  A 

9-4  a  2 

a  -(r13XV13} 

1/2 

13T1 

”r13 " V13~ 

-4  a 

2  -4  a  2 

0  -(r23xv23) 

1/2 

23T2 

rs  —IT  •  V  + 

23  v23 

Equation  (5.18a)  requires  that  the  first  collision  has 
occurred,  and  that  it  is  a  single  overlap  collision  (i.e. 
r23<a)  *  T^e  (5.18b)  require  that  the  second 


(5.19) 


50 


AEDC-TR -72-142 


collision  occurs  at  some  after  the  first  collision. 

The  time  12  is  between  the  (12) -collision  and 

the  separating  (23) -collision.  Condition  (5.18c)  requires 
that  particles  1  and  3  collide  before  2  and  3  separate, 
i.e.  the  (13) -collision  is  also  a  single  overlap  collision. 

In  the  original  integral  (5.13)  the  integration  extends 
over  the  entire  phase  space  of  the  three  particles.  The 
6-function  in  the  T-operator  guaranteed  that  the  integrand 
would  vanish  for  those  phase  points  for  which  the  required 
collision  sequence  would  not  occur.  However,  in  (5.16) 
we  have  integrated  over  these  6-functions,  so  that  the 
integration  must  now  be  restricted  explicitly  to  those 
phase  points  for  which  the  required  collision  sequence 
will  occur. 

Next,  we  consider  the  (SN)-  and  (NS) -terms  in  (5.12), 
which  are  represented  by  the  diagrams  in  Fig.  7b, c.  The 
contributions  of  the  (SN)-  and  (NS) -diagrams  are  related 
through  the  symmetry  property  (4.19),  which  implies 


^ sn  ^  NS 


(5.20) 


The  expressions  for  the  matrix  elements  are 

I^/XIsn  +^'^NS 

(SN) 

+4o.  X<Vv2)*Ao,  *(vl'v3>} 


(5.21) 
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The  volume  element  dfJ  is  given  by  (5.17)  and  the  conditions 
for  the  integration  region  now  refer  only  to  the  SN- 
collision  sequence  as  determined  by  the  diagram  in  Fig.  7b, 
or  explicitly,  by  the  conditions  (5.18a,b)  and 

t1  >  t2  (5.22) 

The  diagrams  in  Fig.  7a, b  show  clearly  that  the  conditions 
(5.18a,b)  apply  to  both  diagrams,  while  the  condition, 

smaller  or  larger  than  t 2,  distinguishes  whether  particles 
1  and  3  will  collide  before  or  after  particles  2  and  3 
separate. 

The  corresponding  matrix  elements  [ij;,x32  for  the 
self  diffusion  can  be  obtained  from  (5.16)  and  (5.21) 

by  replacing  tMv^Vj)  and  xCv^v^)  by  iM^)  and  x^), 
respectively. 

The  collision  integrals  for  y=3  in  (5.10)  can  be 
expressed  as  a  sum  of  three  terms 

[^,X3  33)  =  ^,X]R3)  +  [^rXlc3)  +  ^^x3^3)  ,  (5.23) 


where 


"h*]R3)  =  J 

> 

dv1dx2dx3$  (123)1^  (v1,v2)  *T12G°T13G°T12X  (^1'^2^ 

(5.24a) 

-  J 

dv1dx2dx3$(123)ij^  (v1,v2)  *T12G0T^3G0T32x  (v3 ,v2)f 

(5.24b) 

W.x)^31  -  - 

j'  dv1dx2dx3$  (123)  ^  (v^  ,v2)  * 

T12G°T13G°T13Go,r32X  ^V3  'V2*  * 

(5.24c) 
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The  subscripts  refer  to  the  recollision,  cyclic  collision 
and  hypothetical  collision  sequences  as  represented  by  the 
diagrams  in  Fig.  8. 

After  reduction  the  recollision  term  becomes 


[^,X]r3)  =  f  dfl$(123)A0  ^(v1,v2)*Aa*  x(v£,v£) 

(R)  12 

where  the  volume  element  is  given  by 


(5.25) 


d «  =  °2|^12*°12ld°12°2^13 

A 

’°13 

'  da13dTdvldv2dv3 

(5.26) 

The  integration  region  for 

the 

recollision  is  determined 

by  the  diagram  in  Fig.  8a, 

or. 

explicitly,  by  the  conditions 

?12-“l2  >  0  ’  ^13-S13 

<  0 

r  T  >  0  ; 

(5.27a) 

^i2*2*’!^i2^2-*ovi.2*2  <  0 

■ 

t 

tr  >  0 

(5.27b) 

where  we  have  defined 

°*2  -  ?12(2)+?i2TR  ! 

r12 

(2)  =  312+?12t  ; 

-*•  .-*■  A  -  A 

V1  *  vl-  v13*°13  °13  ! 

=  v2  ' 

(5.28) 

=  -r. 


r12 v12* 


Equation  (5.27a)  states  the  condition  for  the  first  two 
collisions  to  occur,  where  the  integrating  variable  t  is 
the  time  between  the  first  and  second  collision.  Equation 
(5.27b)  gives  the  requirements  for  the  last  collision  to 
occur,  where  t_  is  the  time  between  the  second  and  third 
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collision.  Primed  velocity  variables  refer  to  the  velocities 
between  the  second  and  third  collision,  and  r^  (2)  is  the 
relative  distance  between  1  and  2  at  the  time  of  the  second 
collision. 

We  proceed  with  the  cyclic  collision  term 


[<P»X]i3)  =  f  dfl$  (123)  A  iMv.  ,v?)*A  X(vl,v')  , 

ii  12  32 

where  the  volume  element  is  again  given  by  (5.26).  The 
integration  region  for  the  cyclic  collision  is  determined 
by  the  diagram  in  Fig.  8b,  or  explicitly  by  the  conditions 
(5.27a)  and 


(5.29) 


^r32 v32^  ”^ov32  <  0 


Tc  >  0 


(5.30) 


The  quantities  used  in  (5.29)  and  (5.30)  are  defined  in 
terms  of  integrating  variables  by 


?32  = 

?32(2,+?32tC 

r32(2)  a12"a13+V12T 

7 

. 

^  ^  A  A 

V3  = 

v3+(v13-°13,013  ! 

< 

to¬ 

ll 

< 

to 

• 

■9 

\  r  2 

r  A  2 11/2 

32TC 

=  -r32(2)-v-2-[o 

-{  r32  (2)x  v32}  J 

(5.31) 


The  last  term  in  (5.23),  represented  by  the  hypothetical 
collision  diagram  in  Fig.  8c,  reads 


[♦,Xl 


(3) 

H 


=  -f 

(H) 


dft$ (123) A 


'12 


♦  (v1,V2)*A£ 


32 


X(v. 


->■ 

#V, 


(5.32) 


The  volume  element  is  again  given  by  (5.26)  and  the  integration 
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region  for  the  hypothetical  collision  is  determined  by  the 
conditions  expressed  by  the  diagram  in  Fig.  8c,  or  explicitly 
by  (5.27a)  and 

{?32(2)xv32}2-(av32)2  <0  ;  th  >  0  .  (5.33a) 

tr  >  t±  .  (5.33b) 

The  above  quantities,  expressed  in  integrating  variables, 
are 


a32 

v32tH 


,  (2)+v32th 

;  r32  (2)  "  al2_a13+V12T 

• 

9 

■>  A 

2  ~  2l^2 

•?32(2)-v32- 

La  -fr32(2)xv32)  J 

• 

# 

v13Ti  "  ”2a13*v13 


(5.34) 


Note  that  (5.33b)  expresses  the  condition  that  the  particles 
1  and  3  separate  before  the  (23) -collision,  where  is  the 
time  particle  3  spends  inside  the  action  sphere  of  particle 
1.  If  we  would  reverse  condition  (5.33b) ,  i.e.  take 
Tjj  <  we  would  recover  the  diagram  in  Fig.  7c  and  its 

corresponding  (NS) -matrix  element. 

The  corresponding  contributions  for  the  self  diffusion 

are 


-/  dat(l23> 

(R)  12 


i=l 


'12 


»X3 c"3 ^  =  f  dIi${123)Aa  ^(v2)*Aa  x(v^) 

(C)  12 


'32 


(5.35) 
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[*,X]i3>  =  -  j  dfi$(  123)  A  1|>{V  )*A  x(v2) 

H  (H)  °12  2  °32  2 

In  order  to  obtain  the  collision  integrals  involving 

four  successive  collisions  we  substitute  (5.11)  for  y=4 

into  (5.10)  and  use  (4.16).  We  thus  obtain 

t^,X]]3>  =  I^X]rc)  +  ^^](!r>  +  ^/X]^>  +  [»P,x]^r>  ,  (5.36) 

where 


"  J  dy1dx2dx3*(123)«(v1,v2)*T1260Tj;3G0Tj;2G0T32x'(v3rv2)  r 

[^X](13)  =  J  dv1dx2dx3<Ml23)iMv1,v2)*T12G0T5;3G0T^2G0T13x(virV3)  f 

t^,xl^3)  =  -  f  dv1dx.2dx3Q(l23)il)(vllv2)* 

T12G°T13GoT12GoT12GoT23X  ^v2 ,v3* •  (5.37) 

[^'X^hr^  =  ~J  dv1dx2dx3$(123)ip(v1,v2)* 

T^2^o ^i3Go ^i 3G° T 32G0 Ti 3 X rV3)  • 

The  corresponding  diagrams  are  given  in  Fig.  9.  The  contri¬ 
butions  of  the  CR-diagram  and  the  HR-diagram  are  related 
to  those  of  the  RC-diagram  and  RH-diagram  through  the 
symmetry  property  (4.21) 


^'X^CR*  = 

[*,X]<3)  =  (X/^]^) 


(5.38) 
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After  integrating  out  the  6-functions  in  (5.37)  and  using 
(5.38),  the  matrix  element  for  the  RC-collision  reduces 
to 


-  J  «0*tl23)  k,  *<v1.v2).i  Xtv^.vj)  + 
(RC)  L  12  32 


+\2X(51^2,*40321'(?3'72)]  ’ 


32 


(5.39) 


The  volume  element  dfi  is  again  given  by  (5.26).  The 
integration  region  for  the  (RC) -collision  sequence  is 
determined  by  the  diagram  in  Fig.  9a,  or  explicitly  by  the 
recollision  conditions  (5.27)  and 


{?32(3),^2}2-(ov"2)2  <  0 


*RC  >  0 


(5.40) 


The  above  quantities  can  be  expressed  in  terms  of  integrating 
variables  through  (5.28)  and  in  addition 


°32  r32 (2)+v32tRC 

;  r32^3^  “ 

f32(2)+^32TR 

?32U)  =  312-^13+?12t 

i 

v2  =  v2+  v12  *C,12  °12 

'  v3  =  v3 

/ 

v32TRC  =  -r32 <3) *v32- 

'o2-{  f32(3)^52} 

2y2 

(5.41) 


The  RH-collision  sequence  gives  the  following  matrix 
element 
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'*'X1rh)  +  ^'X'hR>  =  -/«1»U23)[a  *(51.«2)*4  a 

/dtj\  12  J2 


(RH) 


+V  ,x(VV*V  * 


'12  32 

with  dfl  as  given  by  (5.26).  The  integration  region  for 
the  RH-collision  is  determined  by  the  diagram  in  Fig.  9c, 
or  explicitly  by  the  recollision  conditions  (5.27)  and 


(5.42) 


{J32(3)kv^2}2-(ov^2)2  <  0 


T  >0 
RH 


(5.43) 


In  order  to  define  all  symbols  completely  in  terms  of 
integrating  variables,  we  need  in  addition  to  (5.28)  the 
relations 


a32  "  r32 (3)+V32TRH  ? 


r32(3)  "  r32 (2)+V32TR 


r32  (2)  =  al2-al3+V12T  ; 


V32TRH  r32*3)’V32_ 


9  9lV2 

a  -{r32  (3)x  vj2) 


(5.44) 


Note  that  we  have  not  explicitly  imposed  the  condition 
>  t|,  where  v^2x.![  «=  as  was  ^one  tlle 

comparable  case  (5.33b) .  The  reason  is  that  contributions 
from  the  RH-region  where  t_„  <  x  automatically  vanish, 
since  we  have  shown  that  sequences  of  four  successive 
collisions  cannot  contain  any  overlap  collisions  (c.f. 
remark  following  (4.11)). 

The  corresponding  contributions  for  the  self  diffusion 
are  obtained  from  (5.39)  and  (5.42)  by  replacing  ty(v^,Vj) 
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and  x(v.,v.)  by  ^(v9)  and  x(vo)  respectively. 

We  remark  that  the  recollisions,  cyclic  collisions  and 
hypothetical  collisions,  defined  by  the  diagrams  of  Fig.  8 
have  a  close  resemblance  to,  but  are  nevertheless  different 
from,  similarly  named  collision  sequences  introduced  in 
earlier  papers  [3,8].  In  the  present  report  we  have 
separated  the  statistical  and  dynamical  correlations,  so 
that  the  C-  and  H-  collisions  in  Fig.  8b  and  8c  do  not 
contain  any  overlap  collisions,  in  contrast  to  the  collision 
sequences  derived  originally.  Furthermore,  in  the  earlier 
formulation  the  diagrams  representing  collision  sequences 
had  to  be  supplemented  with  a  number  of  auxiliary  conditions 
[3,8].  Those  conditions  effectively  accounted  for  the 
contributions  from  sequences  of  four  collisions.  In  the 
present  report  we  have  shown  that  the  latter  contributions 
can  be  reduced  to  two  collision  integrals  defined  in  terms 
of  the  diagrams  Fig.  9a  and  9c. 
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VI  Discussion 

We  have  shown  explicitly  how  the  first  density  corrections 
to  the  transport  coefficients  of  a  hard  sphere  gas  are 
determined  by  a  set  of  three-particle  collision  integrals. 

These  collision  integrals  account  for  the  correlations  in 
the  position  and  velocity  variables  of  three  molecules. 

The  correlations  are  of  both  a  statistical  and  a  dynamical 
nature.  The  statistical  correlations  refer  to  correlations 
in  configuration  space  irrespective  of  the  velocities  of  the 
particles.  For  a  gas  of  hard  spheres  they  reduce  to  excluded 
volume  effects  and  lead  to  the  consideration  of  overlap 
collisions.  In  an  overlap  collision  the  center  of  at  least 
one  of  the  two  colliding  particles  is  constrained  to  lie 
inside  the  interaction  sphere  of  the  third  particle.  These 
statistical  correlations  are  represented  mathematically  by 
Mayer  functions  (and  they  are  the  only  correlations  to  be 
considered  in  the  density  dependence  of  the  equilibrium 
properties.)  In  addition  we  need  to  consider  dynamical 
correlations  brought  about  by  sequences  of  successive  collisions. 
These  correlations  involve  both  the  position  and  velocity 
variables  and  can  be  represented  mathematically  by  sequences 
of  binary  collision  operators . 

We  have  made  a  systematic  analysis  of  the  statistical 
correlations  versus  the  dynamical  correlations.  This 
analysis  led  to  a  decomposition  of  the  collision  integrals 
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in  terms  of  integrals  associated  with  the  dynamics  of 
one ,  two ,  three  and  four  successive  binary  collisions 
among  three  particles.  In  order  to  accomplish  this,  the 
triple  collision  operator  was  represented  by  a  binary  collision 
expansion.  The  terms  in  the  binary  collision  expansion  were 
then  ordered  such  that  successive  terms  in  the  series  corre¬ 
spond  to  a  decreasing  number  of  Mayer  functions  and  increasing 
number  of  binary  collision  operators^.  The  resulting 
set  of  collision  integrals  is  summarized  in  Table  III  for 
easy  reference. 

(3) 

As  the  first  term  in  this  ordering  we  recover 

the  coefficient  predicted  by  the  theory  of  Enskog.  This 
collision  integral  involves  the  dynamics  of  only  one  collision, 
while  the  positions  are  correlated  by  the  requirement  that 
both  colliding  particles  must  overlap  with  the  third  particle. 
This  term,  therefore,  reduces  to  the  familiar  binary  collision 
integral  associated  with  the  Boltzmann  equation,  multiplied 
with  the  volume  common  to  the  action  spheres  of  two  colliding 
particles . 

The  second  term  accounts  for  the  effect  of  dynamical 
correlations  brought  about  by  sequences  of  two  successive 
collisions.  In  addition,  the  positions  are  now  correlated  by 
the  requirement  that  only  one  of  the  two  colliding  particles 
overlaps  with  the  third  particle.  This  term,  therefore, 
accounts  for  a  combination  of  excluded  volume  effects  and 

^Note  that  the  sequence  TaG0Ta  counts  as  only  one  collision, 
since  the  operators  refer  to  the  penetrating  and  separating 
part  of  the  same  non-interacting  collision. 
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TABLE  III 

Summary  of  three-particle  collision  integrals 


V 

(number  of 
successive 
collisions) 

contributing 
collision  integrals 

equation 

diagram 

1 

(5.7) 

(Enskog  theory) 

2 

l4'#X]gS^ 

(5.16) 

Fig.  7a. 

(5.21) 

Fig.  7b. 

3 

['KxlJ'35 

(5.25) 

Fig.  8a. 

[^.x]^35 

(5.29) 

Fig.  8b. 

[4-.x1h3) 

(5.32) 

Fig .  8c . 

4 

[*.XlS)  +  l*.X]g) 

(5.39) 

Fig.  9a. 

i*»x1bh)+i*'X1hr 

(5.42) 

Fig.  9c. 
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dynamical  effects.  The  corresponding  collision  sequences 
are  represented  by  the  diagrams  in  Fig.  7,  to  which  we  have 
referred  as  SS-,  SN-  and  NS-collisions .  In  the  SS-collision , 
represented  by  Fig.  7a,  both  collisions  are  overlap  collisions 
and  in  the  SN-  and  NS  collisions,  represented  by  Fig.  7b 
and  7c,  only  one  of  the  two  collisions  is  an  overlap  collision. 

Correlations  caused  by  sequences  of  three  successive 
collisions  are  accounted  for  by  collision  integrals  associated 
with  recollisions,  cyclic  collisions  and  hypothetical 
collisions,  defined  by  the  diagrams  of  Fig.  8.  These  collision 
sequences  do  not  contain  any  overlap  collisions. 

Finally,  the  effect  of  correlations  due  to  sequences  of 
four  successive  collisions  is  determined  by  collision 
integrals  defined  with  respect  to  the  diagrams  of  Figs.  9a 
and  9c . 

We  remark  that  the  same  three-particle  collision 
integrals  were  derived  earlier  in  AEDC-TR-71-51  using  geometrical 
considerations  [11,22].  That  derivation  was  based  on  the 
cluster  expansion  of  the  pair  distribution  function  developed 
by  Green  and  Cohen  [23,24].  An  expansion  of  the  three- 
particle  collision  integrals  in  terms  of  the  number  of 
successive  collisions  was  then  obtained  by  decomposing  and 
rearranging  the  triple  collision  integrals  obtained  earlier 
using  a  surface-integral  method  [3,25].  However,  that 
derivation  requires  rather  intricate  geometrical  considerations 
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and  we  prefer  the  more  direct  algebraic  derivation  presented 
in  this  report.  Moreover,  an  attempt  was  made  to  give  a 
complete  self  contained  presentation,  so  that  the  results 
do  not  depend  on  too  many  previous  papers  and  reports. 

An  attempt  to  separate  the  statistical  and  dynamical 
correlations  was  made  recently  by  Henline  and  Condiff  [9] . 

As  the  first  term  they  also  recover  the  Enskog  contribution. 

In  order  to  assess  the  contribution  from  single-overlap 
collisions,  they  introduce  a  term  called  EVD.  This  abbrevi¬ 
ation  indicates  that  the  term  incorporates  both  excluded 
volume  and  dynamical  effects.  However,  in  deriving  this 
EVD  term,  Henline  and  Condiff  consider  only  collision  sequences 
in  which  the  particles  overlap  initially.  As  a  result 
their  EVD  term  does  not  account  for  all  collision  sequences 
that  involve  single-overlap  collisions.  It  corresponds  to 
our  collision  integrals  [i|i  ,x3  gg^  +^'Xlg^  >  but  does  not 
include  The  latter  term  is  equal  to  IXfWg^ 

as  shown  in  (5.20) ;  it  is  considered  by  Henline  and  Condiff 
implicitly  in  a  term  called  TCI.  The  NS-collision 
sequence,  represented  by  Fig.  7c,  is  obtained  from  the  SN- 

collision  sequences  by  time  reversal,  so  that  the  collision 

(3)  (3) 

integrals  [^»XlgN  and  t^/Xl^g  have  the  same  physical 
origin.  The  matrices  of  collision  integrals  associated 
with  the  EVD  and  TCI  terms  of  Henline  and  Condiff  are  not 
symmetric.  In  this  report  we  have  shown  how  a  systematic 
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analysis  leads  to  a  decomposition  of  the  three-particle 
collision  integrals  into  matrices  that  are  symmetric. 

The  relationship  between  our  collision  integrals  and 
those  of  Henline  and  Condiff  is  shown  in  Table  IV.  Henline 
and  Condiff  did  not  consider  collision  integrals  accounting 
for  four  collisions  involving  three  molecules. 

We  have  analyzed  and  calculated  the  various  collision 
integrals,  derived  in  this  report,  for  all  three  transport 
coefficients.  A  detailed  account  will  be  presented  in  a 
subsequent  technical  report. 
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TABLE  IV 

Comparison  with  the  collision  integrals 
of  Henline  and  Condiff. 


Henline  and  Condi££ 

This  work 

EVD 

TCI 

<H"x]<3)  +  [*,x]<3) 

TC2 

['J'-x)c3> 

TC3 

[*,xli3) 

68 


AEDC-TR-72-142 


Appendix  A.  Collisional  transfer  contributions. 

In  the  main  text  of  this  report  we  have  demonstrated 

KK  KK 

how  the  kinetic  contributions  X  ,  n  and  D,  up  to  terms 
linear  in  the  density  n,  could  be  obtained  from  a  cluster 
expansion  of  the  time  correlation  functions  (2.1).  In  this 
appendix  we  show  how  the  potential  energy  contributions 
X  +  X  and  n  +  r\  can  be  calculated  with  the  same 
technique . 

The  potential  energy  contributions  (2.8)  to  the 
currents  J  contain  derivatives  of  the  intermolecular  potential. 
For  hard  spheres  it  is  convenient  to  eliminate  these  deriva¬ 
tives  by  writing  the  current  J(t)  as  the  time  derivative 
of  a  moment  M(t)  [29] 

• 

J(t)  =  M(t)  (A .  1) 


with 


i=l 

The  term  U(r. .)r.  in  (A. 2)  vanishes  for  a  gas  of  hard  spheres. 
iD  1 

Relation  (A.l)  enables  us  to  relate  the  current  to  pseudo¬ 
streaming  operators.  First,  we  note  that  the  correlation 
function  <J(t^)J(t2)>  depends  only  on  the  time  difference 
(t^tj)  '  SO 
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<J(t1)J{t2)>  =  <M(t1)M(t2)>  =  -<M(t1)M(t2)> 


(A. 3) 


and,  secondly,  we  have  in  accordance  with  (2.15a) 


WM 


(t)  =  WLjM(t)  «  LNWLNM(t). 


(A. 4) 


The  correlation  function  can  now  be  expressed  as 

Jvr 

(A. 5) 

By  means  of  (2.10)  and  (A.l)  we  can  introduce  pseudo  currents 


<J 


Zjj  N  tL 

TdxN  ]]  ^(vi)MlNW(12*  *  *N)e  NlNM* 
N  '  J  i-1 


LNM  =  jk  +  jU 


=-tM  _  tk  _u 
-LnM  -  J  -  J. 


(A. 6) 


with  the  kinetic  parts  given  by  (2.7),  and  the  oollisional 
transfer  parts 

N  N 


i  =  tTT  T.  .im(v?-v2)r.  .  -)  hU£. 

+ ,  X  4  2  l  y  13  f—i  1 

i^j  i=l 

N 


(A. 7) 


'  i£j  13  ••  3  13  *s 


Expressions  for  can  be  obtained  from  (A. 7)  by  replacing 
T^  by  TC\jT  1  from  table  (1.8). 

Again  we  can  separate  the  time  correlation  functions 
in  KK-,  KU-,  UK-  and  UU-parts .  The  KU-  and  UK-parts  are 
equal,  as  can  be  verified  with  the  help  of  the  relations 
in  Table  I .  In  the  present  theory  the  UU-part  of  the  time 
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correlation  function  can  be  neglected,  since  it  contributes 
2 

only  in  0(n  )  to  the  transport  coefficients.  Hence,  to  the 

relevant  order  in  the  density  we  may  write  down  a  new  time 

KK  KU 

correlation  function  for  A=A  +2A  and  Section  II  applies 
immediately . 

Instead  of  (2.17)  we  have  expressions  for  the  complete 
A  and  r]  (to  the  relevant  order  in  n)  which  can  be  obtained 
from  (2.17)  by  replacing  the  functions  ^  by  new  functions 
satisfying 

T' (v,e)  =  T  (v,e)  J(v)+a(v,e)  ,  (A. 8) 


with 


na(v1,e)  =  2  I dx^1  II  *  (v±)U(12  . .  .£)  . 

£= 1  1-2 

00  N  (A.  9) 

•I  TiTinT /dxH'i  n  *<vj)W(12...N)JU  , 

N=£  j  =  £+l 

where  JU  is  given  by  (A. 7) . 

To  lowest  order  in  the  density  only  the  term  (£=1) 
contributes  to  a(v;e).  It  gives  for  the  heat  conductivity 


-* 


e) 


2G0  -|-hUv1+n  /dx2^(v2)W(12)T12  J.(vJ-vJ)?12  j  - 
G0|-TTa3nJx  (v^  ,  (A. 10) 


and  for  the  viscosity 
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an*vl'e*  =  2G0 


n  I 

=  G0j|Tra3nSn(v1)  . 


dx2(t>{v2)W{12)T12 


2"fUV12rl2  } 


{A. 11) 


In  the  derivation  of  the  integrals  (A. 10)  and  (A. 11)  we 

have  used  (2.6),  (2.7),  (2.12),  Table  (1.1)  and  hU  =  pU/n  = 

2  3  2  U 

=  jiro  nkT  +  0(n  )  where  p  =  nkT  +  p  is  the  hard  sphere 

pressure.  In  principle,  the  term  (2,=2)  in  (A. 9)  could  also 

contribute  to  0(n),  but  its  contribution  vanishes  due 

to  (1.2). 


Next,  we  substitute  (A. 10)  and  (A. 11)  in  (A. 8),  multiply 
(A. 8)  on  the  left  with  r  ^ ,  and  obtain  integral  equations, 

— v  ^  ^  ^ 

similar  to  (2.27)  in  which  J^(v)  is  replaced  by  (I+^tto  n)J^(v) 
and  (v)  is  replaced  by  (l+~-ira3n)  ^  (v)  .  The  final  results 
for  the  transport  coefficients  are  then 


,  , .  ,4  3  .  •>  .  ,KK 

A  =  (1+^ttcj  n)A0+nA^ 

8  3  K 

n  =  (1+j-^ira  nJno+nn-L 


(A. 12) 


The  terms  A0  and  n0  are  the  Boltzmann  values  for  the  hard 

sphere  transport  coefficients,  the  terms  proportional  to 

nA0  and  nn0  are  the  collisional  transfer  contributions, 

and  are  identical  to  those  obtained  from  Enskog's  theory. 
KK  KK 

The  terms  nA-j^  and  nn1  are  the  contributions  from  triple 
collisions,  calculated  in  the  main  text. 
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Appendix  B.  Reduction  of  triple  collision  integrals. 

In  this  appendix  we  give  the  details  of  the  reduction 
of  the  triple  collision  integrals  to  the  explicit  forms 
listed  in  Section  V.  As  prototypes  of  this  reduction 
we  treat  two  examples,  namely  the  matrix  elements  corre¬ 
sponding  to  the  SS-collision  and  to  the  recollision. 

We  start  with  the  matrix  element  (5.13a),  which  can 
be  written  as 


=  -  /a?1ax2ax3t(i23){T+2*(v:L,v2)}: 


f23G°f23T13X ^Vl/V3^  ' 


(B.l) 


where  we  have  taken  the  hermitian  adjoint  of  the  T-operator 

and  used  the  fact  that  Tt  commutes  with  0(123)  due  to  (1.7). 
According  to  (2.12),  (1.8)  and  (5.15)  one  can  write. 


=  °2./  dS12l712-S12l\2'<'(5l'72)S(J12-312l 

?12-Si2  >  0 

TUX(?1,?3)  =  o2J  do13  I  v13  •  S13 14  x  (v2  ,v3)  a  (r13-o13)  . 


,  S'  I  S\ 


(B.2) 


^13*°13  <  0 


We  further  recall  that  G0  =  (e-L„) 


13 


-1 


00 

-J 


die  eTeTL° ,  and  that 


the  limit  e+0  is  taken  in  (B.l).  Inserting  the  above 
expressions  in  (B.l)  and  using  (2.5)  we  find 
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[\P»X], 


=  -I  dv1dv2dv3d?12d?13$(123 )c2f  d512|v12*ai2 


*12-S12  '  0 


5  (r12-o12)  A  J  (vx  ,v2)  *0  (o-r23)  fdzQ  (a-  |?23+v23t  | )  . 
12  0 

2  /'dCT13lV13‘a13l6(r13+V13T_013)Aa,,X(Vl,V3) 


Vi3-5i3  <  0 


(B .  3 ) 


In  obtaining  (B.3),  we  use  the  fact  that,  for  any  function 


F (r . . ,v . . ) , 
'  13  ' 


e  f (rij+v1jT'vij>  • 

The  appearance  of  6 -functions  allows  us  to  carry  out  some 
integrations.  Using  the  first  6-functions  in  (B.3)  we 
integrate  over  r^2,  so  that  r^2  has  to  be  replaced  every¬ 
where  by  °12.  The  second  6-function  Will  be  used  to  carry 

/s 

out  the  integrations  over  T  and  a13.  In  order  to  do  so  we 
use  the  relation 


(B.4) 


/  dxcr2  f  da. . |v. . *a. . | 6 ^  (r . .+v. .x-o. .)F(a. . ,x)  = 

0J  J  13  1  13  ID1  13  13  13  13 

00 

/dTa2/daij|?iJ-a.j|6<1)  (t-Tu)  i  .xCij)P(?ij+5i.T,t) 


v.  .  *a. .  <  0 
13  13 


=  S(d-|ri:1KOij|)0tTv)P(oij,T1J) 


(B.5) 


where  and  x^  on  the  last  line  of  (B.5)  are  explicitly 
given  by 
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c. .  =  r . . +v . . t 
ID  ID  ID  V 


-*  -  [  2  -  .  2l1//2 

vi jTy  =  ~rij‘vij~|_a  J 


(B.6) 


In  the  first  step  in  (B.5)  the  3-dimensional  6-function 
is  decomposed  into  a  1-dimensional  6-function  for  the 
component  of  r^j  parallel  to  and  a  2-dimensional  6- 

->•  4- 

function  for  the  components  of  r^ ^  perpendicular  to  . 

In  the  second  step  the  integrals  over  the  time  x  and  the 
solid  angle  a.  .  are  carried  out  to  yield,  respectively,  the 
conditions  0(x^)  and  0  (o-  |  r^  *  v^  ^  |  )  ;  for  the  a  .^-integration 
we  use  the  variable  transformation  £  =  to  obtain 


°2  L 


A  ,  — S\  ,  (  O  1  —f  -f-  A 

da .  .  v.  . ‘o .  .  6  '  {r.  .*v.  .-a.  .* v.  .)  = 

ID  1  ID  ID  1  '  ID  ID  ID  ID 


ID  1  ID 

v.  .  *a.  .  <  0 
ID  ID 


13  13  13  13 


(B.  7) 


-I 


d(2)£6(2)  (S-J..XV..)  =  eta-lr.-xv.-l)  . 
|b[  <  a 


After  carrying  out  the  above  integrations  we  obtain  for  the 
relevant  matrix  element 


(3) 

SS 


- -f  avx 


dv2dv3a2|v12-512[dS12dr13 


*0(v12*a12)0{a-r23)0(o-|r23+v23T1|)0(T1)0(a-|r12xv12| ) . 


•  t  (1231  (vx  ,v2  )  .A^X  (vx ,v3) 


(B .  8) 


where 
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13  “  ?13+?13T1 

7  r23  “  r13“a12 

7 

->  A 

(B.9) 

13xl  =  -r 13 " V13~ 

_c  -(r13xv13)  J 

These  results  can  be  put  in  the  form  (5 . 16) - (5 . 19}  of  the  main 

|  <  a  is  equivalent 
with  <  ^2 •  where  is  defined  in  (5.19).  The  reduction 
of  the  SN  -term  proceeds  along  similar  lines. 

As  a  second  example  we  treat  the  recollision  term, 
defined  in  (5.24a).  We  start  again  by  taking  the  hermitian 
adjoint  of  T^  an d  use  relations  similar  to  (B.2) .  The 
recollision  term  can  then  be  written  as 


•  •  »  -► 
text,  by  realizing  that  |r 


12'  23 


'x1R3>  =  '/dVW?12d?134  (123)°2/d512l?12‘°12lS<?12_°12)- 


V12*°12  >  0 


00 

!?1  '*2>  *  I  ix°2  f  dS13  I  ?13  •  °13 ' 5  <J13+? 


13T'013>- 


v13  *°13  <  0 


00 

•  I ds°2  f  d012  ^12*  d12 I S (i12+?12T+?i2s"S12> i0*2x(?i'?2> 


?'2-S*2  <  0 


(B.10) 


where  the  prime  on  the  velocity  variables  represents  the 

action  of  the  operator  R  ,  where  v! .  =  vf -v!  and 

*13  ID  i  3 


=  R  v. 

°13  1 

ll 

<  + 

(v13*al3)a13  ' 

-  “► 

-¥ 

v' 

=  R  v_ 

—  v_ 

2 

a13  2 

2 

(B . 11) 
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The  first  two  6-functions  are  used  to  carry  out  the  integra¬ 
tions  over  r^2  and  ^3*  we  introduce  the  notation 


rl2(2)  =  r12+v12i  ai2+vi2T 


(B.12) 


the  last  line  in  (B.10)  has  precisely  the  form  of  (B.5) ,  so 

A 

that  we  can  carry  out  the  integrals  over  s  and  a*2>  The 
result  is 


W,X]r3) 


=  f  dv1dv5dv^a2  jv10,a-15|da19a2|v7-1‘d1^[dd,-,dT 


23  1  12  w12'  u12u  1  13  13 1  13 

•  0 (vl2 • a12 ) 0 (x) 0 (-v13 *  o13  >  0 (tr) 0 (a- 1 r12 (2) x v|2 | )  . 


(B.13) 


•  $  (123 )  A  ip(v,  fv,)*A  *  iMv',v'), 
a12  x  z  a12  1  * 


where  all  quantities  appearing  in  (B.13)  are  expressed  in 

terms  of  integrating  variables  in  (5.28).  Equation  (B.13) 

,  \ 

is  equivalent  with  the  results  (5.25) - (5.28)  of  the  main  text. 

The  remaining  collision  integrals,  corresponding  to 
the  C-,  H-,  RC-  and  RH-collision  can  be  reduced  in  a 
similar  way. 
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